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Abstract 

Numerous kinds of uncertainties may affect an economy, e.g. economic, political, and 

environmental ones. We model the aggregate impact by the uncertainties on an econ- 
omy and its associated financial market by randomised mixtures of Levy processes. We 
assume that market participants observe the randomised mixtures only through best 
estimates based on noisy market information. The concept of incomplete informa- 
tion introduces an element of stochastic filtering theory in constructing what we term 
"filtered Esscher martingales". We make use of this family of martingales to develop 
pricing kernel models. Examples of bond price models are examined, and we show 
that the choice of the random mixture has a significant effect on the model dynamics 
and the types of movements observed in the associated jdeld curves. Parameter sensi- 
tivity is analysed and option price processes are derived. We extend the class of pricing 
kernel models by considering a weighted heat kernel approach, and develop models 
driven by mixtures of Markov processes. 

Keywords: Pricing kernel, asset pricing, interest rate modelling, jdeld curve, ran- 
domised mixtures. Levy processes, Esscher martingales, weighted heat kernel, Markov 
processes. 

1 Introduction 

In this paper, we develop interest rate models that offer consistent dynamics in the short, 
medium, and long term. Often interest rate models have valid dynamics in the short term, 
that is to say, over days or perhaps a few weeks. Such models may be appropriate for the 
pricing of securities with short time-to-maturity. For financial assets with long-term maturi- 
ties, one requires interest rate models with plausible long-term dynamics, which retain their 
validity over years. Thus the question arises as to how one can create interest rate models 
which are sensitive to market changes over both short and long time intervals, so that they 
remain useful for the pricing of securities of various tenors. Ideally, one would have at 
one's disposal interest rate models that allow for consistent pricing of financial instruments 
expiring within a range of a few minutes up to years, and if necessary over decades. One 
can imagine an investor holding a portfolio of securities maturing over various periods of 
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time, perhaps spanning several years. Another situation requiring interest rate models that 
are valid over short and long terms, is where illiquid long-term fixed-income assets need 
to be replicated with (rolled-over) liquid shorter-term derivatives. Here it is central that 
the underlying interest rate model possesses consistent dynamics over all periods of time 
in order to avoid substantial hedging inaccuracy. Insurance companies, or pension funds, 
holding liabilities over decades might have no other means but to invest in shorter-term 
derivatives, possibly with maturities of months or a few years, in order to secure enough 
collateral for their long-term liabilities reserves. Furthermore, such hedges might in turn 
need second-order liquid short-term protection, and so forth. Applying different interest 
rate models validated for the various investment periods, which frequently do not guar- 
antee price and hedging consistency, seems undesirable. Instead, we propose a family of 
pricing kernel models which may generate interest rate dynamics sufficiently flexible to 
allow for diverse behaviour over short, medium and long periods of time. 

We imagine economies, and their associated financial markets, that are exposed to a 
variety of uncertainties, such as economic, social, political, environmental, or demographic 
ones. We model the degree of impact of these underlying factors on an economy (and finan- 
cial markets) at each point in time by combinations of continuous-time stochastic processes 
of different probability laws. When designing interest rate models that are sensitive to the 
states an economy may take, subject to its response to the underlying uncertainty factors, 
one may wonder a) how many stochastic factor processes ought to be considered, and b) 
what is the combination, or mixture, of factor processes determining the dynamics of an 
economy and its associated financial market. It is plausible to assume that the number of 
stochastic factors and their combined impact on a financial market continuously changes 
over time, and thus that any interest rate model designed in such a set-up is by nature time- 
inhomogeneous. The recipe used to construct interest-rate models within the framework 
proposed in this paper can be summarised as follows: 

(i) Assume that the response of a financial market to uncertainty is modelled by a family 
of stochastic processes, e.g. Markov processes. 

(ii) Consider a mixture of such stochastic processes as the basic driver of the resulting 
interest rate models. 

(iii) In order to expUcitly design interest rate models, apply a method for the modelling 
of the pricing kernel associated with the economy, which underlies the considered 
financial market. 

(iv) Derive the interest rate dynamics directly from the pricing kernel models, or, if more 
convenient, deduce the interest rate model from the bond price process associated 
with the constructed pricing kernel. 

The set of stochastic processes chosen to model an economy's response to uncertainty, 
the particular mixture of those, and the pricing kernel model jointly characterize the dy- 
namics of the derived interest rate model. We welcome these degrees of freedom, for any 
one of them may abate the shortcoming (or may amplify the virtues) of another. For exam- 
ple, one might be constrained to choose Levy processes to model the impact of uncertainty 
on markets. The fact that Levy processes are time-homogeneous processes with indepen- 
dent increments, might be seen as a disadvantage for modelling interest rates for long 
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time spans. However, a time-dependent pricing kernel function may later introduce time- 
inhomogeneity in the resulting interest rate model. The choice of a certain set of stochastic 
processes implicitly determines a particular joint law of the modelled market response to 
the uncertainty sources. Although the resulting multivariate law may not coincide well with 
the law of the combined uncertainty impact, the fact that we can directly model a particular 
mixture of stochastic processes provides the desirable degree of freedom in order to con- 
trol the dynamical law of the market's response to uncertainty. In this paper, we consider 
"randomised mixing functions" for the construction of multivariate interest rate models 
with distinct response patterns to short-, medium-, and long-term uncertainties. Having 
a randomised mixing function enables us to introduce the concept of "partially-observable 
mixtures" of stochastic processes. We take the view that market agents cannot fully ob- 
serve the actual combination of processes underlying the market. Instead they form best 
estimates of the randomised mixture given the information they possess; these estimates 
are continuously updated as time elapses. This feature introduces a feedback effect in the 
constructed pricing models. 

The reason why we prefer to propose pricing kernel models in order to generate the 
dynamics of interest rates, as opposed to modelling the interest rates directly, is that the 
modelling of the pricing kernel offers an integrated approach to equilibrium asset pricing in 
general (see Cochrane [13J, Duffie [15j), including risk management and thus the quan- 
tification of risk involved in an investment. The pricing kernel includes the quantified total 
response to the uncertainties affecting an economy or, in other words, the risk premium 
asked by an investor as an incentive for investing in risky assets. In this work we first 
consider a particular family of pricing kernel models, namely the Flesaker-Hughston class 
(see Flesaker & Hughston UTTl, Hunt & Kennedy [22J, Cairns [llj, Brigo & Mercurio [6]). 
Since our goal in this paper is to primarily introduce a framework capable of addressing 
issues arising in interest rate modelling over short to long term time intervals, we apply 
our ideas first to the Flesaker-Hughston class of pricing kernels. We conclude the paper by 
introducing randomised weighted heat kernel models, along the lines of Akahori et al. [TJ 
and Akahori & Macrina []2l], which extend the class of pricing kernels developed in the first 
part of this paper. 



2 Randomised Esscher martingales 

We begin by introducing the mathematical tools that we shall use to construct pricing ker- 
nel models based on randomised mixtures of Levy processes. We fix a probability space 
(£7, ^, P) where P denotes the real probability measure. 

Definition 2.1. Let {Lf}f>o be an n-dimensional Levy process with independent components, 
and let X : D. ^ M"" be an independent, m-dimensional vector of random variables. For 
t,u e M_,_, the process {M^^(X)} is defined by 

exp{Ku,X)L,) 

^'"^^^-E[exp(h(u,X)LO|X]' ^^-^ 
where thefixnction h : M+ x M'" ^ M" is chosen such that E [ | M^^QC) | ] < oo for all t^ R+. 
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Proposition 2.1. Let the filtration {^t}t>o given by = a ({Ls}o<s<t>^)- Then the 
process {Mf„(X)} is an ({J^f.} ,¥)-martingale. 

We note that X is ^o'^i^^surable and therefore, that {^[} is an initial enlargement 
of the natural filtration of {Lf } by the random variable X. Furthermore, Mq^[X) = 1 and 
Mc^(X) > for all t, u e M+. 

Proof. The condition that E [ | Mfu(X) | ] be finite for all < t < oo is ensured by definition. 
It remains to be shown that 

E[M,,m\-^s]=M,um (2.2) 



for all < 5 < t < oo. We observe that the denominator in ( [2.1[ ) is ^o'^^^surable so that 
we can write 

^ E\exp(h{u,X)Lt)\J^,] 
E iM.JX) I = r \ \ n • (2.3) 

Next we expand the right-hand-side of the above equation to obtain 
E [exp [h[u,XXLt - Lj] exp [h(u,X)Lj | J^f,] 



E [exp [h(u,X)(Lt - Lj] exp [h{u,X)L,'\ |X] 



(2.4) 



Given X, the expectation in the denominator factorizes since Lf — is independent of L^. 
In addition, the factor exp[h(u,X)L5] is J^fJ -measurable so that we may write 

L'^"^^'^^ E[exp[/T(u,X)L,] IX] E[exp[h(u,X)(Lf-Lj] IX] ■ ^ ^ ^ 

Since the increment — and X are independent of L^, the ^-conditional expectation 
reduces to an expectation conditional onX. Thus, equation ( |2.5| ) simplifies to 

exp rh(u,X)Lj 

E [Mf JX) J^fJ = ^ — , Vi 1 ' (2-6) 

L tuv sj E [exp [/i(u,X)Lj] |X] 

which is Msjj(X). □ 

We call the family of processes {Mfy(X)} parameterised by u e the "randomised 
Esscher martingales" (see Gerber & Shiu II20II and Yao [26] for details on the Esscher trans- 
form). The randomization is produced by /i(u,X) which we call the "random mixer". 

Example 2.1. Let {W(-}f>o be a standard Brownian motion that is independent of X, and 



set Lf = Wf in Definition [23] Then, 

Mf JX) = exp [Ku,X)W, - I h\u,X)t] . (2.7) 

Example 2.2. Let {Yt}t>o be a gamma process with rate parameter m > and scale param- 
eter k: > 0. Then E[yf] = Kmt and Var[)'f] = K^mt. We assume that {jt} is independent 



of X. Set Lf = y-f in Definition 2.1 Then, if h(u,X) <k ^, we have 



Mf,(X)= [l-fch(iz,X)]'"'exp[h(iz,X)rf]. (2.8) 
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3 Filtered Esscher martingales 



In this section we construct a projection of the randomised Esscher martingales that can 
be interpreted as follows. Let us suppose that the exact combination of Levy processes 
that forms the stochastic basis of the martingale family {M^^(X)} is unknown. That is, we 
may have little knowledge about how much each of the Levy processes involved actually 
contributes to the stochastic evolution of {M(.y(X)}. The random vector h(u,X) however, 
can naturally be interpreted as the quantity inside {M^^j^X)} that determines at time u the 
random mixture of Levy processes driving the martingale family. Given a certain set of in- 
formation, the actual mixture might not be fully observable, though. This leads us to the 
following construction that applies the theory of stochastic filtering. For simplicity, we focus 
on the case where X is a one-dimensional random variable. 



We introduce a standard Brownian motion {B Jf>o on (CI, ^, P), and define the filtration 

Wt} by 

% = a ({BJo<.<t, {LJo<s<t.^) , (3.1) 



where {Bj is taken to be independent of X and {Lj. We consider the pair 

dXf = 0, 

dIt=tlit,X)dt + dBt, 



(3.2) 
(3.3) 



where £ : M_,_ x M — > M is a well-defined function. The solution to the signal equation ( 3.2 1 



is of course the random variable X. In the theory of stochastic filtering, the process {/t}t>o 
is the so-called observation process. We have 



(3.4) 



It= l{_s,X)ds+Bt. 

Next, we introduce the filtration {^t}t>o defined by 

— ({4}o<s<t> {^s}o<s<t) > 



(3.5) 



where c The filtration {J^J provides full information about the Levy process {L^}, 
however it only gives partial information about the random variable X. Let us thus consider 
the stochastic filtering problem defined by 

M,„ = E[M,JX)|J?t]. (3.6) 

We emphasize thatX is not -measurable and thus {M(.„(X)} is not adapted to {^t}. The 



filtering problem (3.61 is solved in closed form by introducing 

l{_s,X)dB,-l l^{_s,X)ds\, 
Jo Jo J 

where for all t > 



rt 



E 



£{_s,xfds 



< oo. 



(3.7) 



(3.8) 
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and 



E 



r t 



< oo. 



(3.9) 



The process {<?(.} is a ({^^(.},P) -martingale (see, e.g., Bain & Crisan [4]), and it may be used 
to define a change-of-measure density martingale from P to a new measure B by setting 



(3.10) 



The B-measure is characterised by the fact that {/J is an ({^J,B)-Brownian motion. The 
Kallianpur-Striebel formula then states that 



E[M,,(X)|if,] = 
This can be simplified to obtain: 

E[M,JX)|J?J = 



E^ 



(3.11) 



(3.12) 



where the -measurable conditional density /f (x) of the random variable X is given by 

/o(x) exp ( £(5, x)d7, - \ l\s, x)ds') 



(3.13) 



J!l/o(y) exp (Jj£(5,y)d7, - i j^lKs,y)ds) dy' 
A similar filtering system is considered in a different context in Filipovic et al. [H6||. Further 



conditions are imposed on the dynamics of the information process defined in (3.2) and 



(3.3 1, which may be regarded necessary from a modelling point of view. 



Proposition 3.1. Let {Jff } be given hy ^3.5} , and define theprojection = E [M^^iX) \ ^^], 
where {Mfy(X)} is given by (2.11 Then, for t, u e R_^_, {Mf^l is an {{^t},F)-martingale fam- 
ily- 

Proof Recall that c for all t > 0. For 5 < t, we have 

E [Mt, \^s] = E [E [M, JX) I i?,] I , 
= E[M,JX)|J?,], 
= E[E[M,„(X)|^^,] I J?], 



= M„ 



(3.14) 



where we make use of the tower property of the conditional expectation, and the fact that 
{Mtu(X)} is a {^^J-martingale — since c and {Bj is independent of X and {Lj. □ 
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Filtered Brownian martingales. We consider Example |2.1[ in which the total impact of 
uncertainties is modelled by a Brownian motion {Wj. The corresponding filtered Esscher 
martingale is 



(3.15) 



Mtu =\ ftM exp [Ku, x)Wt - lh\u, x)t)dx, 



where the density process {/t(x)}, given in (3.13), is driven by the information process 
defined by p31 ). 

Proposition 3.2. The filtered Brownian models have dynamics 



where 



dM,^ = MMf,{x) [h{u, x)dW, + V,{x)dZ,] dx, 

J —oo 

Mf Jx) = exp [h{u,x)W, - lh\u,x)t] , 
V,{x) = l{t,x)-E[l{t,X)\^,], 



r t 



E [i{s,X)\S^,] ds, 



and /t(x) is defined in (3.13 ). 



Proof. We first show that 

dMtSx) = h{u,x)M,Sx)dW,. 
In Filipovic et al. [|16ll it is proved that 

d/t(x) = (£(t, x) - E [£(t,X) I J ) dZ„ 
where {Z(.}f>o is an ({J2^f},P)-Brownian motion, defined by 

Z, = I,- { E [l{s,X)\^,] ds. 
Jo 

Thus by the Ito product rule, we get 

d[M,Jx)/,(x)] = f,{x)dMM + MMdf,{x) 
since dWf dZf = 0. This simplifies to 



(3.16) 

(3.17) 
(3.18) 

(3.19) 



(3.20) 
(3.21) 

(3.22) 
(3.23) 



d[MMftM]=MMf,ix) 
and we obtain 



f 00 



h{u,x)dW,+ £(t,x)- 



^(t,y)/t(y)dy 



]dz, 



dM, 



where we define 



MMftM [Ku,x)dW, + Vt(x)dZ,] dx 

X) 

/• 00 

v,(x) = £(t,x)- £(t,y)A(y)dy. 



(3.24) 
(3.25) 

(3.26) 

□ 
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Remark 1. The dynamics of {M^^} can be written in the following form: 

dM,^ = E [M,JX)/i(ii,X) I dW, + E [M,JX)V,(X) | JfJ dZ,. (3.27) 

Filtered gamma martingales. Let us suppose that the total impact of uncertainties on an 
economy is modelled by a gamma process {y^} with density 

lP(rt e dy) = ^ \ dy, (3.28) 

where m and k are the rate and the scale parameter, respectively. The associated ran- 
domised Esscher martingale is given in Example 2.2 where h(u,X) < k~^. The correspond- 
ing filtered process takes the form 



f oo 



A(x)([l-;c/i(u,x)]'"'exp[h(iz,x)r,])dx (3.29) 



for h{u, x) <K ^, and where the density /t(x) is given by (3.13 1. 

Filtered compound Poisson and gamma martingales. We now construct a model based 
on two independent Levy processes: a gamma process (as defined previously) and a com- 
pound Poisson process. The idea here is to use the infinite activity gamma process to 
represent small frequently-occurring jumps, and to use the compound Poisson process to 
model jumps, which are potentially much larger in magnitude, and may occur sporadically. 
Let {Cf }t>o denote a compound Poisson process given by 



Ct = Y,^i (3.30) 



i=l 



where {iVt}t>o is a Poisson process with rate A. The independent and identically distributed 
random variables Y; are independent of {iVj. The moment generating function is given by 

E [exp {q Q)] = exp [At (My(e) - l)] (3.31) 

where My is the moment generating function of Yj. For < k"^, we have 

exp(/ii(u,X)n+h2(",W) 



E [exp (hi(u,X)r, + h2{u,X)C,) \X] 

exp exp (h2(i/,X)Cf) 



E [exp {h^iu,X)r,) \X] E [exp Qi2(u,X)C,) \X] 

= mII\x) M^^\XI (3.32) 

where, conditional onX, exp (hi(u,X)7f) and exp {h2(.u,X)Ct) are independent. Further- 
more, 

mII\x) = (l-Khi(u,X))'"'exp(hi(u,X)r,), (3.33) 
M^,^\X) = exp[h2(",X)Q-At(My(h2(",X))-l)]. (3.34) 
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Then, the filtered process takes the form 

r 00 
J —oo 

X exp x) ft + /i2(",^)Q - At (My(h2(",^)) - 1) ] dx, (3.35) 
where ftix) is given by ( |3.13| ). 

4 Filtered Esscher martingales with Levy information 

Up to this point, we have considered a Brovmian information process given by equation 



(3.31. However, the noise component in the information process may be modelled by a 



Levy process with randomly sized jumps, that is independent of the Levy process {Lj used 
to construct the randomised Esscher martingale. In what follows, we give an example of 
continuously-observed information, which is distorted by gamma-distributed pure noise. 

Example 4.1. Let {7t}t>o be a gamma process with rate and scale parameters m and k, 
respectively. We define the gamma information process by 



(4.1) 



Brody & Friedman [T] consider such an observation process in a similar situation. We 
define the filtration {'S^tl by 



and by 



— ({7s}o<s<t' {^s}o<s<t'-^) J 



(4.2) 



(4.3) 



where {/J is given by ( [4.1[ ). To derive the conditional density of X given Jf^, we first show 

(4.4) 



that {If. } is a Markov process with respect to its own filtration. That is, for a e : 

P[7t <a|4,4^,...,4j =P[7,<a|4] 
for all t > 5 > 5i > . . . > 5n > and for all n > 1. It follows that 



[/f < a|7„7,^,...,7,J = 



It <a 



L L 

■'5' T J • • • J r 



Xrt<a 



S S„_i 



Ts 



(4.5) 



It can be proven that /y^, ...,75 /y^ _ are independent of y^ and y^ (see Brody et al. fl9ll). 
Furthermore, fs^^/Ys, Ys„/Ys„_i independent of X. Thus we have 



[7, <a|7„7,^,...,7,J =P[7,<a|7,]. 



(4.6) 



We assume that the random variable X has a continuous a priori density fo[x). Then the 
conditional density of X, 

ftM=-^F[X<x\It], (4.7) 



9 



is given by 

AW = 



Mx)p{l,\X = x) 
r^My)piit\x = y) dy 
/o(x)x--^exp[-/,/(Kx)] 



(4.8) 



where we have used the Bayes formula. The filtered Esscher martingale is thus obtained by 

Mt„ = E[Mt„(X)|^t]. (4.9) 

The result is: 

Q , , /o(x)x-'^^exp[-J,/(i?x)] ^ 

^tu = Mtu(x) -T^ z dx. (4.10) 

5 Pricing kernel models 

The absence of arbitrage in a financial market is ensured by the existence of a pricing kernel 

{Tif }(.>o satisfying > almost surely for all t > 0. We consider, in general, an incomplete 
market and let {St}t>o denote the price process of a non-dividend paying asset. The price 
of such an asset at time t < T is given by the following pricing formula: 

St = — E[7TrSr|JJ-t] . (5.1) 

The price of a discount bond system with price process {PtT}o<t<T <oo ^i^d payoff Pff — 1 
is given by 

1 r 

PtT^—^[nT\^t]- (5.2) 

The specification of a model for the pricing kernel is equivalent to choosing a model for 
the discount bond system, and thus also for the term structure of interest rates, and the 
excess rate of return. A sufficient condition for positive interest rates is that {n^} be an 
({^t},P)-supermartingale. If, in addition, the value of a discount bond should vanish in 
the limit of infinite maturity, then {TZf} must satisfy 

limE[7rr]=0. (5.3) 

T— >oo 

A positive right-continuous supermartingale with this property is called a potential. Let 
{Af}f>o be an {^t}-adapted process with right-continuous non-decreasing paths, where 
Aq — almost surely, and let {AJ be integrable, that is, E [A^] < oo where A^ :— 
limf_>ooAt. Then any right-continuous version of the supermartingale 

Ct = E[A^|if,] -A, (5.4) 
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is a potential of class (D][^ see Meyer ||24l| . Let us denote by {Ct}t>o the potential generated 
by {Af}. Meyer [24] proved that a potential belongs to the class (D) if, and only if, it is 
generated by a process {A^}. Thus, it is enough to choose a process {AJ to model the 
pricing kernel. 

Flesaker & Hughston II17II provide a framework for constructing positive interest rate 
models, in which the pricing kernel is modelled by 



(5.5) 



7rt = J p(u)mtudi;, 

where {nifu}o<t<u<oo is a family of positive unit-initialized martingales, and 

p(t) = -5,Pot. 

It can be shown that the pricing kernel ( [5.5^ is a potential generated by 



r t 



p(u) m^^, du, 



(5.6) 



(5.7) 



and thus, that it is a potential of class (D). Furthermore, given a potential (5.4) where {A^} 
is an increasing, integrable process of the form 



r t 



ds, 



with {at}t>o ^ nonnegative process, there exist a deterministic function 

1 



p(u)=— E[aJ, 

TZq 



and a positive martingale 



(5.8) 



(5.9) 



(5.10) 



for each fixed u> t where mg^ = 1, such that the class (D) potential can be written in the 
form 



f oo 



Tit = Tin 



p{u)m^^ du. 



(5.11) 



where tiq is a scaling factor. Thus, the Flesaker-Hughston models are precisely the class 
of pricing kernels that are class (D) potentials where {AJ is increasing, integrable and of 
the form (5.81, see Hunt & Kennedy [22J. Therefore, to model such class (D) potentials, it 
suffices to specify a family of positive martingales. 

In what follows, we construct explicit Flesaker-Hughston models, which are driven by a 
randomised mixture of Levy processes. We develop such a class of pricing kernels by setting 



f oo 



p{u)M,^du 



(5.12) 



right-continuous {^J-adapted stochastic process {XJt>o is said to belong to the class (D) if the random 
variables are uniformly integrable, where t is any finite {^J-stopping time. 
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where the martingale family {Mfu}o<t<u<oo is defined by (3.61 with > and Mg^ = 1. 
Then, the discount bond system is given by 

PtT = % , , ^ • (5.13) 

The associated instantaneous forward rate {rf7-}o<t<r is defined by r^j = —djlnP^j. We 
deduce that 

rtT = roo , , ^ ' (5.14) 
p(u)Mf„du 

and that the short rate of interest {r Jf>o is given by the formula 

= .00 , , ^ , ^ (5-15) 
Jt p(u)Mf,du 

where :=rtf The interest rate is positive by construction. We note here that the pricing 
kernel models proposed in Brody et al. [8J can be recovered by considering a special case 
of the random mixer, namely h(u,X) = h[u). 



6 Pricing kernel models driven by filtered Brownian martingales 



In the case where the filtered martingales driving the pricing kernel are Gaussian processes, 
the dynamics of the discount bond system can be expressed by a diffusion equation of the 
form (6.2 1. Inserting the filtered Brownian martingale family (3.15 ) into (5.13 ), we obtain 
the price process of the discount bond in the Brownian set-up: 



tT 



J^p(u) j^^ftM exp [Ku,x)W, - lh\u,x)t] dxdu 
I^pM I-ooftM exp [Kv,y)W, - lh^liv,y)t] dydv' 



(6.1) 



A similar expression is obtained for the associated interest rate system by plugging (3.151 
into ( |5.15D . 



Proposition 6.1. The dynamical equation of the discount bond process is given by 

= [rt - 0tt(0tr - dtt) - v„(vtr - Vtt)] dt + (e^r - On)dW, + (v,^ - v,,)dZ, (6.2) 



dPtT 



PtT 

where 



dtT '■ — 



'tT 



p{u)E[MMKu,X)\^,] du 

^ p{u)M,^du 
p{u)E[MMV,{X)\^,] du 
p(_u)Mtudu 



(6.3) 



(6.4) 



T=f 
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Proof. First we have 



p(i;)dMf;j du 



(6.5) 



where dM^^ is given by (3.271. Also, 



p(u)Mt„du 



p(u)dMtudu - p(t)Mff dt. 



(6.6) 



We then apply the Ito quotient rule to obtain the dynamics of {Ptrl- We observe that the 
discount bond volatilities are given by 



(6.7) 
(6.8) 



The market price of risk associated with {W^} is A^^-* := — 0ff; the one associated with |Zf} 
is Af^"" := — V(.(.. The product between the bond volatility vector Qf j. = (ri^^-*, O^^"*) and the 
market price of risk vector = (A^^-*, Ap"*) gives us the risk premium associated with an 
investment in the discount bond, that is, 

n,T ■ A, = -e,, {e,T - e,,) - v,, (v,^ - v„) . (6.9) 

□ 



Proposition 6.2. Let {M^^^iX)} be of the class ^.7) , and let {Mt^} in { [5.M| be given by the 
martingale family ( 3.15} . Then the dynamical equation of the forward rate is given by 



where 



and 



dr^j = [Qfj djd^j + Vf-j- djVf-jl dt — djd^j dW^. — djV^j dZj 
J"p(iz)E[M,JX)h(iz,X)|if,] du 



^ p{u)M,^du 

^'^ p(u)E [M,,(X)V,(X) I J^,] du 
p(u)Mt^du 



(6.10) 



(6.11) 



(6.12) 



where V^{X) is defined by (3.181 



Proof. We apply the Ito quotient rule to (5.141 to obtain the forward rate dynamics. We 
make the observations that 



tT "tT 



and that 



djV^T = rtT V^j 



E[M,t{XMT,X)\^,] 
E[M,r(X)Vt(X)|i^,]" 



(6.13) 

(6.14) 

□ 
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In particular, if we set 



Str = e,T-e,„ (6.15) 
AfT = Vfj--Vtt, (6.16) 



then we can express the risk-neutral dynamics of the forward rate by 

dr^j = [Sfj.5j-Sfj. + Afj-^T-Atj-] dt - Sj-Sfj-dWt - Sj-Afj-dZj, (6.17) 
where {Wjf>o and {Zf }f>o are Brownian motions defined by the Girsanov relations 

dWt = dWt + dt, 

dZt = dZf + Af^dt. (6.18) 



The dynamical equation (6.171 has the form of the HJM dynamics for the forward rate 
under the risk-neutral measure, see Heath et al. II21II . 

Example 6.1. As a first illustration, let us now consider the case in which the information 
process is defined by 

It = (TXt + Bt, (6.19) 

where cr is a positive constant. It can be proven that this is a Markov process (see Brody et 
al. ||10 || ). Equation (6.19 1 is a special case of the path-dependent observation process (3.41. 



2.1 



Let {Wf} be a standard Brownian motion that is independent of X. Then from Example 
we have 

JX) = exp [Ku,X)W, - \ ]x\u,X^t\ . (6.20) 

We suppose that the a priori distribution of X is uniform over the interval (a, b), where 
a > and b > 0. We choose to model the random mixer by 

h(u,X) = cexp(-uX) (6.21) 

where c e M. Here X can be interpreted as the random rate of the exponential decay in 
h[u,X). We obtain the following expressions for the bond price 

r!° p(u) f'' exp \uxL + ce-"^Wf - I ( a^x^ + c^e-^"^) t] dx du 

PtT = oo 1 — ' '^6.22) 

p{u) exp [cry/f + ce-"->'Wt - | (cr^y^ + c^e-^"^') t] dy du 

and the associated interest rate 

p(t) j^exp \axIt + ce~^''W^ - \ (cr^x^ c^e"^'^'') t] dx 

= b • (6.23) 

Jf°° P(") exp [cry/f + ce-"-^ - \ (a^y^ + c^e'^"-)') t] dy 

Since the model is constructed from a single Levy process, it is not — strictly speaking — a 
mixture model as described previously. However, it can be viewed as a kind of two-factor 
Brownian model owing to the presence of the observation process {/J. The bond price 
and the associated interest rate are functions of time and the two state variables and I^. 
Thus, it is straightforward to generate simulated sample paths: 
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Figure 1: Sample paths of discount bond with T = 5 and short rate. We use the filtered 
Brownian model with h(u,X) = cexp(— uX) and X ~ U(a,b). We set a = 0, b = 0.1, 
cr = 0.1, c = 0.5 and Pgt = exp(-0.04t)- 

The parameters a and b influence the rate at which exp(— uX) decays, and together 
with c determine the impact of the Brownian motion {Wj on the bond and interest rate 
evolution. When c is close to zero, the impact of {Wj is very small. For sufficiently large 
values of b — a, a or |c|, the numerical integration in the calculation of the pricing kernel 
may fail to converge. For large values of t, we observe that the sample paths of the short 
rate revert to rg. Thus, there is built-in reversion to the initial level of the short rate. 

7 Bond prices driven by filtered gamma martingales 

Let {yj denote a gamma process with E[7f] = Kmt, and Var[)'t] = K^mt. We consider 
a bond price model based on a pricing kernel that is driven by a family of filtered gamma 



martingales given by (3.29). Then, equation (5.13) for the bond price gives the following 
expression: 

r p(u) J" [1 - KKu,x)r' exp [Ku,x)rt] dxdu 

Jf°°p(v) J^/t(y) [l-Jc/i(v,y)]'"'exp[h(v,y)7f]dydv 

We now investigate this bond price model in more detail, and in particular show the effects 
of the various model components on the behaviour of the bond price. 

Example 7.1. Let the information process {/J, driving the conditional density {/t(x)} be 
of the form 

I, = (TtX + B„ (7.2) 

where X is a binary random variable taking the values X = 1 with a priori probability /o(l), 
and X = with probability /o(0). We choose the random mixer 

h(u,X) = cexp[-bu(l-X)], (7.3) 

where c < and b > 0. Then the expression for the filtered gamma martingale simplifies 
to 

M,,=/,(0)exp(ce-''"rt) (l - Kce"^")"' +/,(l)exp (cr,) (1 - kc)'"S (7.4) 
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where 



/o(0) + /o(l)exp(c77,-ia2t) ' /o(0) + /o(l)exp (a/, - ia2t) ' 

(7.5) 

There are a number of degrees of freedom in this model which have a significant impact 
on the behaviour of the trajectories. In what follows, we analyse the degrees of freedom 
one by one. 

A prion probability: When /o(l) = 0, the diffusion {/J plays no role. The sample paths of 
the discount bond and the short rate are driven solely by the pure jump process. The size 
of the jumps decays over time. As /o(l) increases, there is a greater amount of diffusion in 
the sample paths. Furthermore, there is a higher likelihood of obtaining sample paths for 
which the size of the jumps do not decay over time. If /o(l) = 1, then {M^^} is no longer 
u dependent. This yields a stochastic pricing kernel, but flat short rate and deterministic 
discount bond prices, see Figure |2} 




Figure 2: Sample paths for discount bond with T = 5, and associated short rate. We 
use the Brownian-gamma model with h(u,X) = c exp [— bu(l —X)] where X = {0, 1} with 
m = 0.5, K = 0.5, £7 = 0.1, c = -2,h = 0.03 and Pot = exp(-0.04t)- We let (i) /o(l) = 0, 
aO/o(l) = 0.65 and (uO/o(l) = l- 

Information flow rate a: As the information flow rate increases, the investor becomes 
more knowledgeable at an earlier stage about whether the random variable may take the 
value X = or X = 1, see Figure |3} 

Parameters of the gamma process m and k: The rate parameter m controls the rate 
of jump arrivals. The scale parameter k controls the jump size. 

Parameters of the random mixer b and c: The magnitude of c influences the impact 
of the jumps on the interest rate dynamics. When c = 0, the pricing kernel, and thus the 
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short rate of interest, is deterministic. The sign of c affects the direction of the jumps. 
For < c < jc"^, the short rate (discount bond) sample paths have upward (downward) 
jumps. The opposite is true for c < 0. It should be noted that exp (c exp [— —X)] y^), 
and (l — K'c exp [— biz(l — X)])'"'^ behave antagonistically in c. For large t, one term will 
eventually dominate the other. Thus, for both c > and c < 0, the drift of the short rate 
trajectories is initially negative and then becomes positive for large t, see Figure |4) The 
parameter b determines how quickly the jumps are "killed off". Alternatively, b can be 
viewed as the rate of reversion to the initial level of the interest rate. The interest rate 
process approaches the initial rate more rapidly for high values of b. When b = 0, Mf^ is 
no longer u dependent, and we obtain a stochastic pricing kernel, but flat short rate and 
deterministic discount bond prices, see Figure |5] 




Figure 3: Short rate sample paths for the Brownian-gamma model with h{u,X) = 
cexp [-bu(l -X)] and X = {0, 1}. We choose m = 0.5, k = 0.5, /o(l) = 0.8, c = -2, 
b = 0.03 and Pot = exp(-0.04t)- We set (i) cr = 0.005, (ii) cr = 0.4 and (iii) cr = 1.2. 




Figure 4: Short rate sample paths for the Brownian-gamma model with h{u,X) = 
cexp[-bu(l -X)] and X = {0,1}. We set m = 0.5, k = 0.5, /o(l) = 0.5, cr = 0.1, 
b = 0.03 and Pot = exp(-0.04t)- We choose (i) c = -5, (ii) c = and (iii) c = 1.5. 




Figure 5: Short rate sample paths for the Brownian-gamma model with h{u,X) = 
cexp[-bu(l -X)] and X = {0,1}. We let m = 0.5, k = 0.5, /o(l) = 0.5, cr = 0.1, 
c = -2 and Pq^ = exp(-0.04t)- We choose (i) b = 0, (ii) b = 0.005 and (iii) b = 1. 



17 



Compared to Example 6.1 this model is more robust to variation in the values of the 
parameters. An analysis of the sample trajectories suggests that for large t, the short rate 
reverts to the initial level Tq. 



8 Bond prices driven by filtered variance-gamma martingales 

We let {Lf} denote a variance-gamma process. We define the variance-gamma process as a 
time-changed Brow^nian motion w^ith drift (see Carr et al. [|12|n. that is 

L, = er, + j:B^^ (s.i) 

with parameters e M, S > and v > 0. Here {y^} is a gamma process with rate and scale 
parameters m = 1 /v and k = v respectively, and {B^.^ } is a subordinated Brownian motion. 
The randomised Esscher martingale is expressed by 

JX) = exp [h(iz,X)Lf ] (l - 0vh(iz,X) - lj:^vh\u,X)') , (8.2) 
and the associated filtered Esscher martingale is of the form 

f 00 

M,^ = exp [h{u, x)Lf ] (l - Ovhiu, x) - lT,^vh^{u, x)) dx, (8.3) 



where /t(x) may be given for example by (3.13) or a special case thereof, or by (4.81 



depending on the type of information used to filter knowledge about X. This leads to the 
following expression for the discount bond price process: 

^ _ J~p(u) jr^A(x) exp [Ku,x)L,] (l - evKu,x)-li:W{u,x)y'" dxdu 

jrP(^) J°'^A(y)exp[h(v,y)L,] (l - 0v/T(v,y) - |l]W(v,y))'/'' dy dv' 

(8.4) 



We can also obtain an expression for the short rate of interest by substituting (8.3) into 



(5.15 ). We now present another explicit bond pricing model. 



Example 8.1. We assume that X is a random time, and hence a positive random variable 
taking discrete values {x^, . . . , x„} with a priori probabilities {/oCx^), . . . ,/o(x„)}. We sup- 
pose that the information process {/J is independent of {L J, and that it is defined by 

I, = crXt + Bt. (8.5) 

We take the random mixer to be 

h(u,X) = cexp[-b(u-xf\ (8.6) 

where b > and c e M. We see in Figure |6] that the random mixer, and thus the weight 
of the variance-gamma process, increases (in absolute value) until the random time X, and 
decreases (in absolute value) thereafter. 
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10 20 30 40 



Figure 6: Plot of h(u, X;) for Xi = 2, X2 = 5, X3 = 10 and X4 = 20, where b = 0.015 and 
c = 1 (left) and c = -1 (right). 

The associated bond price and interest rate processes have the following sample paths: 



PtT 




Figure 7: Sample paths for a discount bond with T = 10 and the short rate. We use the 
variance-gamma model with h(u,X) = cexp [— b(u — X)^]. We let 6 = —1.5, S = 2 and 
V = 0.25. We set /o(xi) = 0.2, /qCxs) = 0.35, /oCxg) = 0.35, fglix^) = 0.1 and x^ = 2, 
X2 = 5, X3 = 10, X4 = 20. We choose a = 0.1, c = 0.5, b = 0.015 and the initial term 
structure is Pgt = exp(— 0.04t). 

We observe that over time the sample paths of the interest rate process revert to the initial 
level Tq. However, some paths may revert to Tq at a later time than others, depending on 
the realized value of the random variable X. 

9 Chameleon random mixers 

The functional form of the random mixer h(u,X) strongly influences the interest rate dy- 
namics. The choice of h(u,X) also affects the robustness of the model: there are choices 
in which the numerical integration in the calculation of the pricing kernel does not con- 
verge. So far, we have constructed examples based on an exponential-type random mixer. 
However, one may wish to introduce other functional forms for h(u,X) for which we can 
observe different behaviour in the interest rate dynamics, while maintaining robustness. 
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For instance we may consider a random piecewise function of the form 



KU,X) = gi{u)l^u<X} + g2Ml{u>X} (9.1) 

where gj : — > M for j = 1,2. The random mixer now has a "chameleon form": initially 
appearing to be g^, and switching its form to g2 at X = u. This results in the martingale 
{Mfjj}, and the resulting interest rate sample paths, exhibiting different hues over time, 
depending on the choices of gj (j = 1, 2). We can extend this idea further by considering 
(i) multiple gj, or (ii) a multivariate random mixer of the form 

Ku,X, Y^, Y2) = Y,)l^u<x} + g2(", ^2) W}' (9.2) 

where X > 0, Y^ and Y2 independent random variables with associated information 
processes. In this case, the gj are themselves random-valued functions. Here X can be 
regarded as the primary mixer which determines the timing of the regime switch. The 
variables Y; (i = 1,2) can then be interpreted as the secondary mixers determining the 
weights of the Levy processes over two distinct time intervals. 

Example 9.1. We now present what may be called the "Brownian-gamma chameleon 
model". We consider the filtered gamma martingale family (3.291 in the situation where 
the random mixer h(u,X) has the form 

Ku,X) = cj sm{_a^u)l^^<x} + ^2 exp (-a2") '^{u>x} (9.3) 

where Ci,C2 < and a2 > 0. The information process {/J associated withX is taken to 
be of the form 

It = atX + Bt. (9.4) 

We assume that X is a positive discrete random variable taking values {xj, X2, . . . , x„} with 
a priori probabilities /o(^i), i = 1,2, ...,n. That is, the function h(u,X) will switch once 
from sine to exponential behaviour at one of the finitely many random times. Inserting 
(3.29 ), with the specification (9.3 1, in the expression for the bond price ( |5.13 1, we obtain 

CpC") ZI"=i/t(^i) [1 -Kh{_u,Xi)]"'^ exp [h(u,Xi)Y^] du 
p^^ — il ' ^ I: il I: (9.5) 

jrP(^) i;"=i/t(yi) [l-Kh{v,y{)Y^ exp [h{v,yi)Yi] dv' 
where h{u,Xj) is given by ( [9.3[ ) forX = X;, and 

/o(x,)exp[ax,/,-la\^t] 

JAXi)- rr ^ r r 1 2 2.1 ■ ^ 

Since the sine function oscillates periodically within the interval [—1,1], the integrals in 
(9.51 may not necessarily converge to one value. However, at some finite random time 
u=X, the sine behaviour is replaced by an exponential decay; this ensures the integrals in 
the expression for the bond price converge. Such a behaviour may be viewed as a regime 
switch at a random time. In the simulation below, the analysis of the model parameters is 
analogous to the one in Example |7.1[ It is worth emphasizing nevertheless that (i) the a 
priori probabilities /oC^;), i = 1, 2, . . . , n have a direct influence on the length of the time 
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span during which the sine function in the chameleon mixer is activated, (ii) the magnitude 
of aj determines the frequency of the sine wave, while affects the rate at which reversion 
to the initial interest rate (in the simulation below Vq = 4%) occurs, and(iii) the size of z-y 
determines the amplitude of the sine, and it significantly impacts the convergence of the 
numerical integration. We find that reasonable results are obtained for —k~^ < < k~^. 




Figure 8: Sample paths of discount bond with T = 10 and short rate trajectories. 
We use the Brownian-gamma chameleon model with h[u,X) = sin(aii;)l{„<j)f j + 
C2 exp (— a2i/) '^{u>x\- LetX take the values {x^ = 2, X2 = 5, X3 = 10, X4 = 15} with apriori 
probabilities {/o(xi) = 0.2, /o(x2) = 0.35, /oCxg) = 0.35, f^ix^ = 0.1}. We set m = 0.5, 
K = 0.5, cr = 0.1, Ci = 0.2625, C2 = 0.75, = 0.75, 02 = 0.02 and Pot = exp(-0.04t)- 



10 Model-generated yield curves 

The yield curve at any time is defined as the range of yields that investors in sovereign debt 
can expect to receive on investments over various terms to maturity. For a calendar date t 
and a time to maturity t, we let Y^^^^^ be the continuously compounded zero-coupon spot 
rate for time to maturity t, that is, the map t >-> Y^^^^^. We vwite 

Pt,t+r = e^p{-^Y,^t+r). (10.1) 

Typically, the following yield curve movements are observed: (i) parallel shifts of the yield 
curve corresponding to an equal increase in yields across all maturities; (ii) steepening 
(flattening) of the yield curve, that is the difference between the yields for longer-dated 
bonds and shorter-dated bonds widens (narrows), and (iii) changes in the curvature and 
overall shape of the yield curve. The terms "shift", "twist" and "butterfly" are also used to 
describe these yield curve movements. 



As shown in Figure 10 below, the two-factor Brownian-gamma model set-up in Example 
7.1 is indeed too rigid to allow for significant changes in the shape of the yield curve. For 
/o(l) = 1, the yield curve is flat at all times. For < /o(l) < 1, this model can generate 
flat, upward sloping yield curves and in certain cases, slightly inverted yield curves. The 
variance-gamma model (Figure [12] ) and the Brovraian-gamma chameleon model (Figure 



141 show more flexibility, where changes of slope and different yield curve shapes are 
observed. These model may generate flat, upward sloping, inverted and humped yield 
curves. We emphasise that these classes of models are able to capture all three types of 
yield curve movements. 
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Figure 9: Discount bond curves for the Brownian-gamma model. We let X = {0, 1} with 
/o(l) = 0.3. We let m = 2, K = 0.2, a = 0.1, c = -2, b = 0.03, Pot = exp(-0.04t)- 
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Figure 10: Yield curves for the Brownian-gamma model. We letX = {0, 1} with/o(l) = 0.3. 
We let m = 2, K = 0.2, a = 0.1, c = -2, b = 0.03, Pq^ = exp(-0.04t)- 
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Figure 11: Discount bond curves for the variance-gamma model with h[u,X) = 
cexp[-b(u-Xf]. We let 6 = -1.5, S = 2 and v = 0.25. We set /o(xi) = 0.2, 
/o(x2) = 0.35, /o(x3) = 0.35, /o(x4) = 0.1 and Xj = 2, X2 = 5, X3 = 10, X4 = 20. We 
choose a = 0.1, c = 0.5, b = 0.015 and the initial term structure is Pot = exp(— 0.04t). 
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Figure 12: Yield curves for the variance-gamma model where h(u,X) = c exp [— b(u — X)^]. 
We let e = -1.5, S = 2 and V = 0.25. We set /o(xi) = 0.2, /oCxa) = 0.35, /oCxg) = 0.35, 
/q(x4) = 0.1 and = 2, X2 = 5, X3 = 10, X4 = 20. We choose u = 0.1, c = 0.5, b = 0.015 
and the initial term structure is Pq^ = exp(— 0.04t). 
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Figure 13: Discount bond curves for the Brownian-gamma chameleon model. We let X = 
{xi = 2,X2 = 5,X3 = 10, X4 = 20} with /o(xi) = 0.15, /oCxa) = 0.35, /ofe) = 0.35, 
/q(x4) = 0.15. We let m = 0.5, k = 0.5, u = 0.1, q = -0.4375, C2 = -1.25, = 0.75, 
a2 = 0.02, Pot = exp (-0.040- 
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Figure 14: Yield curves for the Brownian-gamma chameleon model. We let X = {x^ = 
2,X2 = 5,X3 = 10, X4 = 20} with /o(xi) = 0.15, /ofe) = 0.35, /ofe) = 0.35, /o(x4) = 
0.15. We let m = 0.5, k = 0.5, a = 0.1, = -0.4375, C2 = -1.25, = 0.75, 02 = 0.02, 
Pot = exp (-0.040- 
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1 1 Pricing of European-style bond options 



Let {Qt}o<s<t<r be the price process of a European call option with maturity t and strike 
< iC < 1, written on a discount bond with price process {Ptr}o<t<r- The price of the 
option at time s is given by 



Cst = —^[n,{P,T-K)+\^s 
By substituting \S.12\ and ( |5.13| ) into ( [ll.ip , we obtain 



(11.1) 



Cst = — E 



1 



C 00 



p(u) du-iC 



p(u)Mt^du 



if. 



(11.2) 



In the single-factor models that we have considered with a Markovian information process 
{/J , we can define the region y by 



r- 00 



C 00 



p(u) = y, /t = z) du - iC 



p(u)MfjLf = y,7t = z) du > 0| . 

(11.3) 



It follows that the price of the call option is 



P(")Mt„(y,z)du-iC p(u)Mt„(y,z)diz q,(y,z)dydz (11.4) 



where 



[L,<y,7,<z|J?,] 



5y 5z 

We can use Fubini's theorem to write this more compactly in the form 



71^ 



p(iz)$tudu-iC 



where 



^tJy,z)qs(y,z)dy dz. 



(11.5) 



(11.6) 



(11.7) 



We apply Monte Carlo techniques to simulate option price surfaces. A large number of it- 
erations is required to obtain accurate estimates. To increase precision, variance reduction 
techniques or quasi-Monte Carlo methods can be considered (see Boyle et al. [5 J). The 
choice of the random mixer affects the shape of the resulting option price surface. The sim- 
ulations in Figure 15 are based on (i) the Brownian-gamma model constructed in Example 



7.1 and (ii) the Brownian-gamma chameleon model in Example 9.1 The wave across 



the second option price surface is produced by the sine function that defines part of the 
chameleon random mixer. 
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Figure 15: Option price surface at s = 2 of call options on a discount bond with T = 10. 
(i) Simulation based on the Brownian-gamma model. We setX = {0, 1} with/o(l) = 0.5, 
m = 0.5, K = 0.5, a = 0.1, c = -2, b = 0.03 and Pgt = exp(-0.04t). (ii) Simulation based 
on the Brownian-gamma chameleon model. We set X = {x^ = 2, X2 = 5, X3 = 10, X4 = 20} 
with/o(xi) = 0.15, /o(x2) = 0.35, /ofe) = 0.35, /o(x4) = 0.15, m = 0.5, k = 0.5, u = 0.1, 
Ci = 0.35, C2 = 1, aj = 3, a2 = 0.03, and Pq^ = exp(— 0.04t). 



12 Randomised heat kernel interest rate models 

In Sections 2 and 3, we constructed martingales based on Levy processes and an Esscher- 
type formulation. We recall that the pricing kernel is modelled by 

= 



t 

r 00 



p(u)E[Mf,(X,Lj li^,] du 

p(u)M,, (x,L,) du/t(x)dx. (12.1) 
The process {M^^ (X, Lf)} is a unit-initialized positive {'^J -martingale, and the process 

r 00 

S, (X,L,) := p(i/)Mf„ (X,Lf) diz (12.2) 



is a positive {'S^J-supermartingale. The projection of a positive {'^f }-supermartingale onto 
{JfJ, that is 

Tif :=E[St(X,Lj IJZ^f], (12.3) 
is an {<^J-supermartingale (Follmer & Protter II18II . Theorem 3). 

Weighted heat kernel approach. We now model the impact of uncertainty on a finan- 
cial market by a process that has the Markov property with respect to its natural filtration, 
and which we denote {Y^}^>q. Of course, the case where {Yf.} is a Levy process, which is a 
Markov process of Feller type, is included (see Applebaum [l3[|). 

Definition 12.1. Let {Y^} be a Markov process with respect to its natural filtration. A mea- 
surable function p : M_,_ X M_,_ X M — > M 15 a propagator if it satisfies 



E[p(t,v,yJ \Y,] =p{s,v + t-s,Y,) 



(12.4) 
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for (v, t) e M+ X and < 5 < t. 



Next, let {"f }f>o be a pure noise process, and let the filtration {^^J be generated by 

= ({^Jo<.<t> {"Jo<.<t, X) , (12.5) 

where {Y^}, {n^}, and the random variable X are all independent. Let G(-) be a positive 
bounded functioE|^ and let h : M+ x M ^ M. Then we set 



p(t, V, := E [G (h(t + v,X), 7^+,) | . 
This is a {'^J -propagator since X is 'S^g'^i^^surable. It follows that 



s,(x,yj:= 



w{t,v)E[Gm + v,XlY,+,)\%] dv 



(12.6) 



(12.7) 



is a {^f}-supermartingale, see Akahori et al. [T\. Here w(t,v) is a positive function that 
satisfies 

w(t,v-5) <w(t-5,v) (12.8) 
for arbitrary t, v e M_,_ and s < t A v. Now we define the market filtration {.^J by 

SO<s<t' [Is SO<s<t ) , (12.9) 

where {/J carries information about X, which is distorted by the pure noise }. We have 
that c '^f. Then, by Follmer & Protter [18J Theorem 3, the projection 



n, ■=E[S,{X,Y,)\^,] 
is an {J^f }-supermartingale. It follows that 



(12.10) 



E 





f 00 



w(t,v)E[G(h(t + v,X),y,+,) 1^,] dv 
w(t,v)E[E[G(h(t + v,X),yt+v)l^t] li^t] dv. 
w(t,v)E[G(h(t + v,X),yt+v)l^t] dv. 



(12.11) 



We emphasize that in equation (12.11 1, E [G(h(t + v,X),Yf_f_^) is not an {^J -propagator 
when {/J is not a Markov process. Nevertheless, {tt,- } is a valid model for the pricing ker- 
nel, subject to regularity conditions. 



^Once a Markov process {Y,} has been chosen, it may be sufficient to relax the boundedness condition, and 
choose G(-) to be a positive and integrable function. 
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13 Quadratic model based on the Ornstein-Uhlenbeck process 

In this section, we generate term structure models by using Markov processes with de- 
pendent increments. We emphasize that such models cannot be constructed based on the 
filtered Esscher martingales. Let us suppose that {Y^} is an Ornstein-Uhlenbeck process 
with dynamics 

dYf = 5(/3 -yJdt-FTdWf, (13.1) 

where 5 is the speed of reversion, (3 is the long-run equilibrium value of the process and T 
is the volatility. Then, for s < t, the conditional mean and conditional variance are given by 



E[yjy,] = exp[-5(t-s)] + ^(l-exp[-5(t-5)]) 



Var[yjy,] = — (l-exp[-25(t-5)]). 
Let us suppose, for a well-defined positive function h : M_,_ x 

G(h(v,x),yj = h(v,x)y2. 



that 



(13.2) 
(13.3) 

(13.4) 



Since X is ^^o'^i^asurable, and by applying (13.2 1 and { 13.3 ), it follows that 

pliu,t,Y„X) = E[h(t + u,X)y2^J^t] , 

= h(t + u,x)E[(n+„-E[y,+jy,]+E[y,+jy,])2 \y, 

= h{t + u,X) [Var in] + E [Y,+^ \Y,f], 

(l - e-2^") + [Yt e"^" + /3 (l - e"^") ] ' 



= h{t + u,X) 



25 



(13.5) 



The pricing kernel is then given by (12.11 ), and we obtain 



Tit = 



w(t,u) 



'Y^Z 

25 



(l - e-25") + [y, e"^" + I3{l- e"^") ] " 



X 



h(^t + u, x)/f (x) dx du. 



(13.6) 



It follows that the price of a discount bond is expressed by 



= — E 



Tit 



w{T,v)E[G{h{T + v,XlYj+,) dv 







(13.7) 



where {tiJ is given in (13.61, and the conditional expectation can be computed to obtain 



w(r,v) 



"Y^z 
25 



(1 



-25{T+v-t) 



) + [y, e-'^^+'-'^ + /3 (l - e-^(^+^-^)) ] ■ 



KT + v,x)f,Mdxdv. (13.8) 
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Example 13.1. We assume that X is a positive random variable that takes discrete values 
{xj, . . . , x„} w^ith a priori probabilities {/o(^i), • • • ,/o(^n)}- We suppose that the information 
flow^ {7f } is governed by 

I, = aXt + Bt. (13.9) 



We choose the random mixer to be 

h{t + u,X) = Ci exp [-C2(t + u -X)] (t + u), 
where > and C2 > 0, and w^e assume that the weight function is 

w(t,u) = exp [-;(u + t)] 



(13.10) 



(13.11) 



for j > 0. Later, in Proposition |14.1[ we show that this model belongs to the Flesaker- 
Hughston class. Therefore, the short rate of interest takes the form 



r. = 



e-^'E[Gm,XlY,)\^,] 



e-^-^^+^) E [cm + v,Xl I JfJ dv 



(13.12) 



Next we simulate the trajectories of the discount bond and the short rate process. We 
refer to lacus II23II for the simulation of the Ornstein-Uhlenbeck process using an Euler 
scheme. We observe oscillations in the sample paths owing to the mean-reversion in the 
Markov process. 




Figure 16: Sample paths for a discount bond with T = 10 and the short rate for the 
quadratic OU-Brownian model with /i(t + u,X) = exp (— C2(t + u — X))(t + u) with 
Ci = 0.02 and C2 = 0.1. We let 5 = 0.02, p = 0.5, T = 0.2 and Yq = 1- We let = 1 and 
X2 = 2 where /o(^i) = 0.3 and /o(^2) = 0-7 and a = 0.1. The weight function is given by 
w(t,u) = exp [-0.04(t + u)]. 



The model-generated yield curves follow. In this example, we mostly observe changes of 
slope and shifts. However, it should be possible to produce changes of shape in the yield 
curve by varying the choices of G(-) and h(-). 
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5 10 15 20 25 30 

Figure 17: Discount bond curves for the quadratic OU-Brownian model with h(t + u,X) = 
Ciexp(-C2(t + u-X))(t + u) with = 0.01 and C2 = 0.1. We let 5 = 0.02, /3 = 0.5, 
T = 0.2 and Yq = 1. We let = 1 and X2 = 2 where /o(xi) = 0.5 and /oCxa) = 0.5 and 
cr = 0.1. The weight function is given by w(t,u) = exp [— 0.04(t + u)]. 




Figure 18: Yield curves for the quadratic OU-Brownian model with h(t + u,X) = 
Ciexp(-C2(t + u-X))(t + u) with = 0.01 and C2 = 0.1. We let 5 = 0.02, /3 = 0.5, 
T = 0.2 and Yq = 1- We let = 1 and X2 = 2 where /o(xi) = 0.5 and /oCxa) = 0.5 and 
cr = 0.1. The weight function is given by w(t,iz) = exp [— 0.04(t + u)]. 
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14 Classification of interest rate models 



In what follows, we show that, under certain conditions, the constructed pricing kernels 
based on weighted heat kernel models belong to the Flesaker-Hughston class. 

Proposition 14.1. Let {Y^} be a Markov process, and let the weight function be given by 
w(t, v) = o/^Ct + v) with i/) : ^ M+ such that 



^{s)E[G{h{s,X),L,)] ds < oo. 



(14.1) 



Then, the pricing kernel 



V^(t + v)E[G (h(t + v,X),Lt+,) \ ^t] dv 



(14.2) 



15 a potential generated by 



^{u)G{h{u,X),L^) du. 



(14.3) 



that is, a potential of class (D). Thus, the pricing kernel is of the Flesaker-Hughston type. 



Proof. The function w(t, v) = i/)(t + v) satisfies ( 12.8 ), and thus is a weight function. Then 
we see that 



n, = 



= ^0 



where 



i/)(iz)E [G(h(iz,X),Lj] m,,du 
pMm^^ du, 

E[G (h(iz,X),Lj IJf,] 



(14.4) 



(14.5) 



E[G(h(iz,X),L,)] 

is a positive unit-initialized {J^j }-martingale for each fixed u > t. The constant tzq is a 
scaling factor. □ 

We note that, for instance, the potential models of Rogers [25] which can be gener- 
ated by the weighted heat kernel approach with xp(t + v) = exp [— a(t -I- v)] where a > 0, 
are Flesaker-Hughston models. To generate potentials from the weighted heat kernel ap- 
proach with a general weight w(t,v), the weight function and G(-) should be chosen so 
that E[7rt] ^ as t ^ oo. 

Let us suppose that {yj is a Markov process with independent increments. Then the 
class of Esscher-type randomised mixture models presented in this paper, for which 



exp [h(u,X)Lf] 
E [exp [h(u,X)Lf] \X] ' 



(14.6) 
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cannot be constructed by using the weighted heat kernel approach. We see this by setting 



G{h{v,X),L,+,) = 



exp [h(v,X)Lt+^] 
E [exp [h(v,X)L,+,] \X]' 



(14.7) 



and by observing that E[G(h(v,X), Lj+y) | '^f] is not a {'S^f} -propagator. As we mentioned 
earlier, the class of models introduced by Brody et al. dSjl is included in the class of Esscher- 
type randomised mixture models. Similarly, models based on kernel functions of the form 
G(h(x), Ff) can produce other Esscher-type models by use of the weighted heat kernel 
approach. The following is a diagrammatic representation of the considered classes of 
positive interest rate models: 



FLESAKER-HUGHSTON MODELS 




WEIGHTED HEAT KERNEL MODELS 



Figure 19: Classes of interest rate models. 



We conclude with the following observations. The pricing kernel models proposed in 
this paper are versatile by construction, and potentially allow for many more investigations. 
For instance, we can think of applications to the modelling of foreign exchange rates where 
two pricing kernel models are selected — perhaps of different types to reflect idiosyncrasies 
of the considered domestic and foreign economies. In this context, it might be of particular 
interest to investigate dependence structures among several pricing kernel models for all 
the foreign economies involved in a polyhedron of FX rates. We expect the mixing function 
h{u,X) to play a central role in the construction of dependence models. Furthermore, a 
recent application by Crisafi of the randomised mixtures models to the pricing of 
inflation-linked securities may be developed further. 
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